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We study a two-particle circular billiard containing two finite-size circular particles that collide elastically with 
the billiard boundary and with each other. Such a two-particle circular billiard provides a clean example of an 
"intermittent" system. This billiard system behaves chaotically, but the time scale on which chaos manifests 
can become arbitrarily long as the sizes of the confined particles become smaller. The finite-time dynamics 
of this system depends on the relative frequencies of (chaotic) particle-particle collisions versus (integrable) 
particle-boundary collisions, and investigating these dynamics is computationally intensive because of the long 
time scales involved. To help improve understanding of such two-particle dynamics, we compare the results 
of diagnostics used to measure chaotic dynamics for a two-particle circular billiard with those computed 
for two types of one-particle circular billiards in which a confined particle undergoes random perturbations. 
Importantly, such one-particle approximations are much less computationally demanding than the original 
two-particle system, and we expect them to yield reasonable estimates of the extent of chaotic behavior in 
the two-particle system when the sizes of confined particles are small. Our computations of recurrence-rate 
coefficients, finite-time Lyapunov exponents, and autocorrelation coefficients support this hypothesis and 
suggest that studying randomly perturbed one-particle billiards has the potential to yield insights into the 
aggregate properties of two-particle billiards, which are difficult to investigate directly without enormous 
computation times (especially when the sizes of the confined particles are small). 



A traditional billiard system consists of a point 
particle confined in some domain (which is usu- 
ally a subset of R 2 ) and colliding perfectly elas- 
tically against the boundary of that domain; 10 i 26 
Such billiards can have chaotic, regular (i.e., inte- 
grable), or mixed dynamics ,L2i2s For example, a 
finite-size circular particle confined within a cir- 
cular boundary is integrable, but two circular par- 
ticles confined in a circular domain yield chaotic 
dynamics, though regular behavior can persist for 
extremely long times 1 ^. 

The chaotic dynamics in two-particle billiards 
appears via the dispersive mechanism 26 as a result 
of particle-particle collisions, whereas particle- 
boundary collisions lead to regular dynamics. 
Consequently, although the long-time dynamics 
is chaotic, the regular transients can become ar- 
bitrarily long as one considers confined particles 
with progressively smaller radii. It is desirable to 
find means to simplify investigations of the statis- 
tical properties arising from long-time transient 
dynamics in two-particle billiards, whose dynam- 
ics are not well understood and which require 
very long computations to simulate. In this pa- 
per, we take a step in this direction by considering 
one-particle billiards with random perturbations 
and comparing diagnostics for measuring aggre- 
gate levels of chaotic dynamics in two-particle 
versus perturbed one-particle billiards. In par- 
ticular, we consider two circular particles con- 
fined in a circular billiard and one-particle cir- 
cular billiards with two types of random pertur- 



bations: one in which random perturbations are 
applied at times determined via a Poisson pro- 
cess and another in which these times are given 
by actual particle-particle collision times from the 
two-particle system. 

The two-particle circular billiard considered in 
the present paper provides a clean example of an 
"intermittent" system, so this example deserves 
many future investigations. 



I. INTRODUCTION 

Systems of hard particles interacting via perfectly elas- 
tic collisions provide paradigm examples for studying 
the foundations of statistical mechanics^ Among the 
primary examples used to study classical and quantum 
chaos in conservative systems are billiard systems ) 10 ^ 26 
which can be implemented experimentally in both clas- 
sical and quantum situations . n i 18 ] 22 

Billiards are one of the most important types of Hamil- 
tonian system^ Typical classical Hamiltonian systems 
are neither fully chaotic nor fully regular (i.e., integrable) 
but instead have "mixed" dynamics. 11 That is, their 
phase space has both regular and chaotic regions. How- 
ever, generic mixed systems are very difficult to analyze, 
so it is important to study Hamiltonian systems with 
simpler but non-generic mixed dynamics that allow more 
thorough analysis /' 12 ! 19 ' 24 ! 25 Billiard systems are among 
the most important systems for such pursuits. 

The most commonly studied type of billiard system 
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consists of a single point particle confined in a closed 
planar region that collides against the boundary of the 
region such that its angle of incidence equals its angle 
of reflection. In recent years, the study of open billiards 
(which contain a hole or leak through which particles 
can escape) has become increasingly prominent 4 -? . How- 
ever, it is much less common — but nevertheless extremely 
interesting — to study few-particle billiards, in which a 
small number of finite-size confined particles collide elas- 
tically both against a billiard boundary and against each 
otherJ^ The balls move freely in between collisions. 

Depending on the geometries of the billiard boundaries 
and the confined particles, it is possible for few-particle 
billiards to exhibit both regular and chaotic features. 
This situation, however, is somewhat different from the 
idea of mixed dynamics mentioned above. Namely, al- 
though the dynamics of few-particle billiards are fully 
chaotic in the infinite-time limit, it is possible to con- 
struct systems such that integrable dynamics last for ar- 
bitrarily long periods of time. In particular, this can 
occur in two-particle billiard systems in which both the 
billiard balls and the billiard table are shaped like cir- 
cles as one considers confined particles with progressively 
smaller radiiJ^ Two-particle billiards provide a clean ex- 
ample of an "intermittent" system 11 . This is an impor- 
tant motivation for studying this example, as intermit- 
tent systems contain "sticky" regions near which typical 
trajectories can get trapped for very long times and can 
be notoriously difficult to study in detail. Accordingly, 
it is useful to examine "simple" (relatively speaking) ex- 
amples of such systems. 

A circular billiard table with a confined finite-size cir- 
cular particle is integrable, and every ball-boundary col- 
lision in a two-particle circular billiard leads to regular 
dynamics. However, particle-particle collisions in this 
two-particle system ultimately lead to chaotic dynamics 
(via the dispersive mechanism) in the infinite-time limit. 
The dynamics of this system arise from a competition 
between the integrable particle-boundary collisions and 
the chaotic particle-particle collisions. The relative fre- 
quency of the latter versus the former becomes smaller 
as the radii of the confined particles become smaller, and 
the integrable transients can therefore last for arbitrar- 
ily long times. Consequently, although the dynamics are 
eventually chaotic if one waits long enough, one must 
examine the transient dynamics to achieve a thorough 
understanding of this system. 

Investigating the transient behavior of two-particle bil- 
liards entails very long computation times, which is par- 
ticularly true when the integrable transients are long. 
Accordingly, we use numerical computations to attempt 
to discern when a perturbed circular one-particle billiard 
can give a reasonable approximation for statistical prop- 
erties of the two-particle billiard. When this is the case, 
which we show happens for small particle radii, one can 
attempt to gain insights into the original two-particle sys- 
tem by studying a perturbed one-particle system. 

The rest of this paper is organized as follows. In Sec- 



tion nn we describe the model of a circular billiard with 
confined circular particles of finite radius. For simplic- 
ity, we assume that both particles are of the same size. 
In Section IIII1 we describe some measures of the sta- 
tistical properties of the system's dynamics. In Section 
HVi we compare these statistical properties for the two- 
particle circular billiard and for two types of perturbed 
one-particle circular billiards. In Section [Vj we summa- 
rize our results. 



II. MODELLING A CIRCULAR BILLIARD SYSTEM 

We study a two-particle billiard that occupies a re- 
gion QcK 2 with a circular boundary dQ in which each 
finite-size, circular particle of mass m and normal mo- 
mentum p n = mv n moves freely until it encounters ei- 
ther the boundary or the other particle. To understand 
this system, we need to start from the simplest case. Ac- 
cordingly, we assume in this paper that m — \v\ = 1 for 
both particles and also that they have the same radius 
(we consider several different radii in our computations). 1 
Without loss of generality, we set the radius of the bil- 
liard table to be 1. 

Consider the dynamics of one of the two particles in 
our system. When it collides against the boundary, its 
momentum changes according to 

p' =p- 2(p,n)n, (1) 

where p is the momentum before the collision, p' is the 
momentum after the collision, and ft is the unit normal to 
the boundary dQ at the collision point. If the particles 
collide against each other instead, then the change in 
momentum of one particle still obeys (TTJ), but now ft is 
the normal to the tangent of the boundary of the other 
particle at the point of collision. 

In a dynamical system, a Poincare section is a lower- 
dimensional subspace of phase space that is transversal 
to the flow. In a two-particle billiard, we can obtain a 
Poincare section separately for each particle. Consider 
the coordinates (s,-,sin(0.j)), where Sj € [0, 27r] is the ar- 
clength coordinate of the boundary at a collision between 
particle i and the boundary, and 6i G [0, 2tt] is the cor- 
responding angle of collision against the boundary— We 
construct a Poincare section by tracking the coordinate 
values of such collisions for each particle while throwing 
away the continuous flow and particle-particle collisions 
in between such particle-boundary collisions. This yields 
a 4-dimensional Poincare section, as two dimensions from 
phase space are eliminated by only considering the places 
of particle-boundary collisions. (We lose one dimension 
from using a Poincare section and a second dimension via 
energy conservation.) 

In Fig. [T] we show 2-dimensional projections of ex- 
amples of Poincare sections for one of the particles in 
the two-particle system. We consider examples in which 
each confined particle has a radius of 0.02 (top panel) or 
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0.014 (bottom panel). We plot the dependence of sin(^i) 
of particle-boundary collisions versus the boundary loca- 
tion (i.e., the arclength Si) of the collision. To construct 
Fig.[TJ we use discrete-time simulations with 15000 time 
steps in which each time step corresponds to 0.3 time 
units of the simulation. When the radii of the confined 
particles are small, as in these two examples, the dynam- 
ics can exhibit regimes of regular-looking behavior even 
for long run times. 

Whether one can discern transient regular dynamics in 
the Poincare sections depends on the initial conditions, 
computation time, and the sizes of the confined parti- 
cles. In particular, regions of regular behavior (which 
appear as incomplete lines and circles) are more evident 
in the top panel then in the bottom panel, although there 
are regions of regular behavior in the latter plot as well. 
Such regular features arise from particle-boundary colli- 
sions, which are angle-preserving. Horizontal segments 
represent consecutive particle-boundary collisions with 
the same angle, and parts of circles represent regions of 
recurrent behavior. Chaotic features arise from particle- 
particle collisions, whose dispersive character results in 
changes to 9i in subsequent particle-boundary collisions, 
leading to divergence of trajectories and the splattering 
of dots in Poincare sections. 

In these simulations, we determine initial conditions as 
follows. For each of the two confined particles, we choose 
initial distances between its center and the center of the 
billiard table uniformly at random from (0,1). As the 
particles have finite radii (i.e., they are not point parti- 
cles), we discard any configuration in which the particles 
overlap with each other or any particle overlaps with the 
billiard boundary. We do this by checking the overlap- 
ping condition and choosing a new random initial posi- 
tion until there is no overlap with other particle or the 
boundary. We choose the initial angle of each particle 
uniformly at random from the interval [0, 2ir). Because 
the magnitude of the speed for each particle is always 
equal to 1, this gives the starting velocity of each parti- 
cle. 

In this paper, we perform computations for many val- 
ues for the radius r = r\ = r-i of the confined parti- 
cles. For a given family of computations, we use the 
same ensemble of initial conditions for each choice of r. 
For a given family, we thus choose initial conditions as 
described in the above paragraph using the largest em- 
ployed value of r. This guarantees that there is never any 
particle-particle or particle-boundary overlap. 



III. QUANTIFYING CHAOTIC DYNAMICS 

Although the Poincare sections in Fig.Q]illustrate regu- 
lar features that result from particle-boundary collisions 
and chaotic dynamics that result from particle-particle 
collisions, it is desirable to quantify the statistical prop- 
erties resulting from the interplay of these two types of 
events as functions of the sizes of the confined particles. 
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FIG. 1. (Color online) Projections of example Poincare sec- 
tions for one particle in a two-particle circular billiard. We 
use the same initial conditions for both panels. (As discussed 
in the main text, we choose initial conditions uniformly at 
random.) In panel (a), the radii of each of the the confined 
particles is 0.02, and we can observe fairly regular behavior 
(horizontal segments from consecutive particle-boundary col- 
lisions and parts of circles from recurrent behavior). In panel 
(b), the radii of each confined particle is 0.014. The behavior 
now appears to be more chaotic, as there are more isolated 
points and fewer regular features. 

As diagnostics to describe chaotic dynamics in our sys- 
tem, we use recurrence-rate coefficients, autocorrelation 
coefficients, and finite-time Lyapunov exponents. 



A. Recurrence Plots and Recurrence Rates 

Recurrence plots (RPs) were introduced as a tool 
to visualize the recurrences in a variable Xi in phase 
spacei 14 ! 21 Recurrence of states is a typical feature of non- 
linear and chaotic systems, and it is traditional to try to 
find them by visualizing high-dimensional phase spaces 
as projections to two-dimensional or three-dimensional 
spaces. RPs take a different approach and yield a visu- 
alization using an N x N matrix, where N denotes the 
number of time steps in a simulation. 
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To construct an RP, we start with the formula^ 



Ri,j = 9(e 



x j\\) 



(2) 



where Xi G R m and i,j G {1, ...,N}. Additionally, e 
is a threshold distance, || • || is a norm, and 0(-) is the 
Heaviside function. In this paper, we always use the L^- 
norm 28 . 

We compute RPs for the temporal evolution of the 
positions of the confined particles. We choose one of the 
two particles in the system and construct its recurrence 
plot. The variable Xi represents the (x, y) position of 
the chosen particle inside the billiard at time i. We are 
interested in how close and how often a particle returns 
to (x, y) coordinates that it has visited previously in its 
trajectory. An RP includes a dot in the location 
if the positions of the particle center at times i and j 
are within a threshold distance e from each other. The 
threshold should be small enough so that the particle can 
be considered to have approached sufficiently close to the 
previously visited locations, but it should also be large 
enough to keep computations reasonably efficient^ 

In this paper, we present results using the thresh- 
old value e = 0.001. (Computations using the value 
of e — 0.01 gave similar qualitative results.) Early 
RP studies suggested that one should choose e to be a 
few percent of the phase-space diameter—, and similar 
comments have been made about problems with circu- 
lar symmetry^ Additionally, some authors have used an 

-dependent threshold but we have elected to use 
a uniform value. The constant value for the threshold 
yields a symmetry in the recurrence plot: the RP always 
has a main diagonal and is symmetric about it. 

In Fig. [5J we show an example RP for one particle in 
a two-particle circular billiard. As with Poincare sec- 
tions, splatters of dots illustrate irregular behavior. The 
lines parallel to the diagonal represent regions of regular 
behavior, in which trajectories visit the same region of 
the phase space at different times. The isolated points 
by themselves don't contain any information about the 
system. However, the isolated points occurring next 
to the lines in this RP can be an indication of chaotic 
dynamics^ which we know occurs in this system. The 
length of diagonal segments is determined from the du- 
ration of the temporal evolution of regular trajectories.— 

An RP can be used to define a recurrence-rate (RR) 
coefficient 14 ' 21 



RR(e) 



N 2 



(3) 



which measures the density of recurrence points in an 
RP and can be used as a measure of complexity in a 
dynamical system. As the discrete time N —> oo, RR 
represents the probability for a state to recur in an e- 
neighbourhood in phase spaced 
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FIG. 2. (Color online) Example recurrence plot with Eu- 
clidean norm for a two-particle circular billiard with two con- 
fined circular particles (each of which has a radius of 0.02). 



B. Autocorrelation 

Autocorrelations measure the tendency for observa- 
tions of the same system made at different time points to 
be related to one another*^ That is, an autocorrelation 
describes a correlation in a time series with respect to its 
own past and future values. A negative autocorrelation 
value indicates that the direction of influence is changing 
as a function of time, whereas a positive autocorrelation 
value can be construed as a tendency for a system to re- 
main in a similar state from one observation to the next. 

For a series x consisting of N observations, the au- 
tocorrelation coefficient Ri calculated between the time 
series and the same series lagged by I time units is given 
by the formula 

1 = — r^ 2 ' 

where (x) is the mean value of x among the N observa- 
tions. 

In this paper, we will compute autocorrelation coeffi- 
cients for time series of the angle of particle-boundary 
collisions. 



C. Lyapunov Exponents 

The largest Lyapunov exponent A measures the rate of 
separation of trajectories in a dynamical system.— It is 
defined using the equation 



|W(t)|| 



.At I 



M0) 1 1 



(5) 



where Sd(t) is the time-i separation of two trajectories 
that start a distance 8d(0) apart. We use a finite-time 
version of the Lyapunov exponent because of the finite- 
time nature of numerical computations. The dynamics of 
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a system are only predictable up to the Lyapunov time, 
which is defined as the time it takes for two neighboring 
trajectories to diverge by a distance equal to e. 

A positive Lyapunov exponent indicates that trajec- 
tories separate from each other exponentially fast. For 
chaos to manifest, there needs to be a positive Lyapunov 
exponent and trajectories also need to mix. A positive 
Lyapunov exponent implies that there is a local instabil- 
ity, and mixing implies that trajectories of individual par- 
ticles get arbitrarily close to each other arbitrarily often 
if the system evolves for a sufficiently long time . 11 ' 15 ' 27 

Our two-particle system behaves chaotically in the 
t — > oo limits, though its most interesting behavior oc- 
curs during the potentially extremely long transients, and 
we use numerical computations to study such dynamics. 
To calculate the largest finite-time Lyapunov exponent, 
we use Benettin's algorithm^, which assumes that small 
perturbations in initial conditions will stretch primarily 
along the most unstable direction in phase space after a 
sufficiently long time. We evolve the billiard system for 
n length-T time steps. (The time step r thus gives the 
unit of time for our numerical simulations.) This yields 
a finite-time Lypaunov exponent of 

where ||5d(fcr)||/||M(0)|| is the stretching factor due to 
an initial perturbation of size ||<M(0)||. After each of 
the n time steps (which we index by k), we evaluate 
log(||M(fer)||/||5d(0)||) and compute A max as a mean over 
these n evaluations. 

For our simulations, we show results for (r = 0.2, 
n = 25000) and (r = 0.3, n = 15000)^ The size of the ini- 
tial perturbation is ||<W(0)|| = 1/5000. In our numerical 
computations, we report values for the largest finite-time 
Lyapunov exponent A max (n) that is a mean over results 
that we obtain using an ensemble of initial conditions 
(200 or 1000 in the examples shown). We choose these 
initial conditions using the procedure that we discussed 
m Section HH 

We calculate finite-time Lyapunov exponents by sep- 
arately evaluating the stretching factor for the x- 
coordinates and y-coordinates and then taking the max- 
imum of the two corresponding exponents [A^ ax (n) and 
A^(n)] calculated using equation ©. 



IV. COMPARISON OF TWO-PARTICLE BILLIARDS 
AND PERTURBED ONE-PARTICLE BILLIARDS 

Because analytical calculations are very difficult for 
two-particle billiards, we use numerical computations to 
determine when statistical properties for the two-particle 
circular billiard can be approximated effectively by those 
for a perturbed one-particle billiard. Perturbed one- 
particle billiards have the potential to be more tractable 



analytically than two-particle billiards, and they can re- 
quire significantly less computational time (and fewer 
computational resources). 

In this paper, we compare our two-particle billiard sys- 
tem with perturbed one-particle billiards using numerical 
simulations to compute the diagnostics (largest finite- 
time Lyapunov exponents, recurrence-rate coefficients, 
and autocorrelation coefficients) that we discussed in Sec- 
tion mil 

We consider two different types of perturbation to one- 
particle circular billiards, where we apply perturbations 
("kicks") to the confined circular particle. We first exam- 
ine a perturbed one-particle billiard with time intervals 
between kicks determined using the exponential proba- 
bility distribution. This corresponds to using a Poisson 
process to determine the number of perturbations that a 
particle experiences in a given amount of time. As ex- 
plained in Section IIVB1 we estimate the parameter for 
the exponential distribution (and hence for the Poisson 
process) from the statistics of the original two-particle 
system. The second type of perturbed one-particle circu- 
lar billiard that we consider is one with random kicks at 
times that are determined directly from particle-particle 
collisions in the original two-particle billiard. 

A. Random Perturbations in a One-Particle Billiard 

We apply random perturbations to a one-particle circu- 
lar billiard at times determined by a Poisson process (see 
the discussion below) or directly from particle-particle 
collision times in the two-particle billiard. We impose the 
random perturbations by determining a velocity angle 
uniformly at random from the interval [0, 27r). We have 
chosen this type of perturbation because we are consider- 
ing small confined particles, so particle-particle collisions 
are not very common, and we are assuming that we have 
no knowledge of the angle at which such a collision oc- 
curs. We have also considered angles chosen from the 
standard Gaussian N(Q, 1) distribution. This approach 
appears to yield similar qualitative results, so we only 
show the results of our computations for the uniform dis- 
tribution. 



B. Poisson Process for Determining Perturbation Times 

Our original system consists of two circular particles of 
the same size confined inside a radius- 1 circular billiard. 
These particles collide elastically against the boundary 
of the billiard and against each other. As we discussed 
in Section [TTl we choose the initial positions and veloc- 
ities (i.e., angles) of the particles uniformly at random 
for the two-particle billiard. For our initial conditions 
in perturbed one-particle systems, we take the set of ini- 
tial conditions for one of those two particles. We compute 
each of our diagnostics for both two-particle billiards and 
one-particle billiards for each initial condition, and we 
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show plots with the means of those diagnostics over all 
initial conditions. (We consider as many as 1000 initial 
conditions for some calculations.) 

When applying kicks in the perturbed one-particle bil- 
liards, we change the particle velocity (i.e., its angle) as 
explained in Section IIV Al Because of the randomness of 
the initial data, we assume that all of the time intervals 
between particle-particle collisions (and hence between 
particle perturbations in the one-particle systems) have 
the same distribution. Each such collision is thus inde- 
pendent from other collisions. Note that we make this 
assumption for simplicity, and it is desirable to relax it. 

We seek a means to estimate kick times for a particle 
in a perturbed one-particle billiard to attempt to obtain 
an approximation for some statistical properties of the 
original two-particle billiard without having to simulate 
the original two-particle system. 

Because of our independence assumption, we use a 
Poisson process to describe the sequence of particle- 
particle collisions in our two-particle system. (One can, 
of course, use more sophisticated distributions based on 
the system's dynamics, and this is an important idea for 
future investigations.) The Poisson process was devel- 
oped to model events that occur by chance and inde- 
pendently from each other while maintaining a constant 
intensity (i.e., the expected number of events per unit 
time is constant) The Poisson process in our exam- 
ple has a rate of where /x gives the mean (con- 
tinuous) amount of time per particle-particle collision. 
As the times of kicks of the confined particle in a per- 
turtbed one-particle billiard should correspond to the 
times for the particle-particle collisions in the original 
two-particle system, the Poisson process gives estimates 
of these perturbation times. In Fig. [3j we show a plot 
of the perturbation times determined by a Poisson pro- 
cess in a one-particle system versus the original times of 
the particle-particle collisions in the two-particle billiard. 
In this example, the radius of each confined particle is 
r = r\ = T2 = 0.008. We fit a linear function to this 
curve using the method of least squares. The high qual- 
ity of our fit suggests that using a Poisson process to 
determine perturbation times is a reasonable approxima- 
tion. 

It follows from the definition of a Poisson process 
that the time intervals between the perturbations are 
given by independent and identically distributed (IID) 
random variables. Because the system behaves ergod- 
ically in the infinite-time limit, such a distribution for 
particle-particle collisions seems plausible. The sequence 
of particle-particle collisions in the two-particle billiard 
and hence of the perturbations in the one-particle bil- 
liard then behaves as a Markov process, so that the fu- 
ture of the process depends only on the present state and 
is independent of the past. However, when considering 
finite-time dynamics, recurrences in the system lead to 
deviations from the Poisson-process approximation to the 
collision times, and this becomes increasingly noticeable 
as the radii of the confined particles become larger. We 



note, however, that computations of escape rates in circu- 
lar billiards with a hole 9 suggest that the above Poisson 
picture is too simplistic. Examining this in more detail 
is a very interesting idea for future work. 
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FIG. 3. (Color online) Comparison of the times of particle- 
particle collisions in a two-particle circular billiard with times 
of perturbations in a perturbed one-particle system deter- 
mined using a Poisson process for the random initial con- 
ditions and confined particles of radius r — r\ = f2 = 0.008. 

To avoid having to simulate the original two-particle 
billiard to obtain the mean number of particle-particle 
collisions (and hence the rate of the Poisson process) 
that we will subsequently use in a perturbed one-particle 
billiard, we do a set of simulations of the two-particle 
billiard first to obtain estimates of this quantity. In 
Fig. [4] we show the dependence of the number of particle- 
particle collisions on the particle radius in the two- 
particle billiard for the smaller particle sizes that we con- 
sidered. We fit a linear curve to this graph using the 
method of least squares, and the fit is reasonably good. 
This suggest that, for a given particle radius, one can pre- 
determine the rate of the Poisson process for and then 
use that parameter value when investigating a perturbed 
one-particle billiard. However, we did not do this when 
comparing our results for two-particle billiards and per- 
turbed one-particular billiards. Instead, to obtain better 
estimates of the Poisson rates, we perform numerical sim- 
ulations of the original two-particle system for each initial 
condition. 



C. Results of the Comparisons 

In this section, we compute diagnostics for circular 
two-particle billiards (which is computationally inten- 
sive) and compare them to the same diagnostics com- 
puted for perturbed one-particle billiards. For it to be 
reasonable to apply random perturbations to the veloc- 
ity variable in one-particle billiards, we need to consider 
long computations and small confined particles. For sim- 
plicity, we assume that both confined (circular) particles 
in the two-particle billiard have the same size, though of 
course it would be interesting to examine more general 
situations. We consider several different particle radii in 
order to examine the effect of particle size on the diag- 
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FIG. 4. (Color online) Mean number of particle-particle col- 
lisions in a circular two-particle billiard as a function of the 
radius of the confined particles (which both have a radius of 
r = n = r 2 = 0.02). 



nostics. As discussed in Section Hill the diagnostics that 
we use to measure chaotic dynamics are the maximial 
finite-time Lyapunov exponent, the recurrence-rate coef- 
ficient, and the autocorrelation coefficient for the angle 
of particle-boundary collisions. 

In Figs. [5HZJ we show the dependence on the radius of 
the confined particles of the diagnostics applied to one 
particle in our three systems (which we label as 'Two- 
particle', 'Poisson', and 'Actual'). As discussed previ- 
ously, we average each diagnostic over a large number 
of initial conditions. In the plots, we also include error 
bars whose length of 1/5 of the standard deviation corre- 
sponding to the calculated mean value (i.e., the error bars 
show 1/10 above and below the mean). The computed 
standard deviations are large because we choose our ini- 
tial conditions uniformly at random and the systems are 
chaotic. 

As the plots illustrate, the diagnostics have values of 
the same order of magnitude for both the two-particle 
billiard and the perturbed one-particle billiards. In some 
cases, these values are sometimes also very close to each 
other quantitatively. In general, the values of the diag- 
nostics appear to become more dissimilar as one consid- 
ers confined particles of larger radii, though we do not 
observe a monotonic dependence. 

From our numerical computations alone, it is diffi- 
cult to discern when our one-particle systems provide 
"good" approximations to the original two-particle bil- 
liard. However, we believe that examining two-particle 
billiards versus associated perturbed one-particle bil- 
liards has the potential to be valuable. Additionally, 
these perturbed one-particle billiards are interesting sys- 
tems to investigate in their own right. 

Interestingly, drawing the perturbation times from a 
Poisson process sometimes yields better results than tak- 
ing perturbation times from the actual times of particle- 
particle collisions in the original two-particle system, 
even though doing this uses less information from the 
original system. This result, together with Fig. El sug- 
gests that it is reasonable to use a Poisson process to 
estimate the number of particle-particle collisions in the 



two-particle billiard in a certain time interval. 

V. CONCLUSIONS AND DISCUSSION 

Circular two-particle billiards provide a clean exam- 
ple of intermittent dynamics, so studying them in detail 
should be very useful for obtaining a better understand- 
ing of intcrmittency. In the present paper, we have taken 
a step in this direction by comparing the properties of 
a circular two-particle billiard with two different circular 
one-particle billiard systems undergoing random pertur- 
bations. We considered random perturbations at times 
determined in two different ways: (1) times taken from 
a Poisson distribution, and (2) times obtained directly 
from the times of particle-particle collisions in the origi- 
nal two-particle system. 

For such approximations to be reasonable, one should 
consider small confined particles, which entails very 
long computation time for the two-particle system, as 
particle-particle collisions occur much less frequently 
than particle-boundary collisions. Consequently, the 
two-particle system exhibits long transients with regu- 
lar behavior, although the system behaves chaotically in 
the t — y oo limit. Accordingly, it can be very helpful 
to find situations in which it is reasonable to compare 
some of the aggregate properties of two-particle billiards 
to properties in perturbed one-particle systems (which 
require much less exhaustive simulations to study). We 
focused in the present paper on numerical simulations, 
but ideally perturbed one-particle billiards will also be 
studied analytically. 

We considered Lyapunov exponents, recurrence-rate 
coefficients, and autocorrelation coefficients and found 
that the two-particle billiard and perturbed one-particle 
billiards have values of the same order of magnitude when 
the confined particles are small. At times, we also found 
quantitative agreement. As these diagnostics are used 
widely to measure the amount of chaos in a system, this 
suggests that the aggregate levels of chaos in the one- 
particle systems are similar to that in the original two- 
particle system. 

Importantly, the one-particle approximations are much 
less computationally demanding than the original two- 
particle system, and we expect that they will also provide 
a much easier setting for analytical investigations, which 
are very difficult for two-particle billiards. We thus pro- 
pose that randomly perturbed one-particle billiards have 
the potential to provide insights on the aggregate (and in- 
termittent) dynamics of associated two-particle billiards. 
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FIG. 5. (Color online) Finite-time Lyapunov exponents for 
a two-particle circular billiard ('Two Particle'), a perturbed 
one-particular circular billiard with perturbation times drawn 
from a Poisson distribution ('Poisson'), and a perturbed one- 
particle circular billiard with perturbation times drawn from 
the times of particle-particle collisions in the two-particle sys- 
tem ('Actual') as a function of particle radius for (identical) 
confined particles. In the top plot, we consider radii rang- 
ing from 0.008 to 0.03. We average over 200 random initial 
conditions (see the discussion in the text for how we choose 
these initial conditions) and simulate for 25000 time steps (of 
0.2 time units each) for each particle radius. In the bottom 
plot, we consider radii ranging from 0.008 to 0.026. We av- 
erage over 1000 random initial conditions and simulate for 
15000 time steps (of 0.3 time units each) for each particle 
radius. In the top plot, observe that our computations for 
the Poisson-time perturbed one-particle billiard gives a good 
approximation to those for the two-particle billiard. The Lya- 
punov exponents from the Actual-time perturbed one-particle 
billiard have the same order of magnitude as the other two 
sets of values. In the bottom plot, observe that our calcula- 
tions for the Actual-time perturbed one-particle billiard gives 
an excellent approximation to those for the two-billiard. The 
Lyapunov exponents from the Poisson-time perturbed one- 
particle system are also very similar. The mean values are 
located in the centers of the error bars, which depict 1/10 of 
a standard deviation above and below the mean. 
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FIG. 6. (Color online) Recurrence-rate coefficients for a two- 
particle circular billiard ('Two Particle'), a perturbed one- 
particular circular billiard with perturbation times drawn 
from a Poisson distribution ('Poisson'), and a perturbed one- 
particle circular billiard with perturbation times drawn from 
the times of particle-particle collisions in the two-particle sys- 
tem ('Actual') as a function of particle radius for confined 
particles with radii ranging from 0.008 to 0.03. We average 
over 200 random initial conditions and simulate for 25000 
time steps (of 0.2 time units each) for each particle radius. In 
the bottom plot, we show results only for small radii to make 
it easier to see the very good agreement between all three 
computations in that regime. For larger radii, the two per- 
turbed one-particle billiards still have similar RR coefficients, 
but the RR coefficient for the two-particle billiard exhibits an 
interesting balloon (depicted in the top plot) before shrinking 
again. Even during the balloon, all three computations give 
values with of the same order of magnitude. The mean values 
are located in the centers of the error bars, which depict 1/10 
of a standard deviation above and below the mean. 
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FIG. 7. (Color online) Autocorrrelation coefficients for the 
angle of particle-boundary collisions for one particle in a two- 
particle circular billiard ('Two Particle'), a perturbed one- 
particular circular billiard with perturbation times drawn 
from a Poisson distribution ('Poisson'), and a perturbed one- 
particle circular billiard with perturbation times drawn from 
the times of particle-particle collisions in the two-particle sys- 
tem ('Actual') as a function of the (identical) radii of the 
confined particles. In the top plot, we consider radii ranging 
from 0.008 to 0.022. For each radius, we average the auto- 
correlation coefficient over 200 initial conditions and simulate 
for 25000 time steps (of 0.2 time units each). In the second 
plot, we average over 1000 initial conditions and simulate for 
15000 time steps (of 0.3 time units each). We observe ex- 
cellent agreement for all three sets of diagnostics in the top 
plot, and we still find good agreement in the bottom plot 
(though the results for the two-particle versus one-particle 
billiards state to become different quantitatively for larger 
radii). The mean values are located in the centers of the er- 
ror bars, which depict 1/10 of a standard deviation above and 
below the mean. 



